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Using a multiphase transport model, we study the production of multistrange baryons from 
the hadronic matter formed in relativistic heavy ion collisions. The mechanism we introduce is 
the strangeness-exchange reactions between antikaons and hyperons. We find that these reactions 
contribute significantly to the production of multistrange baryons in heavy ion collisions at SPS 
energies, which has been found to be appreciably enhanced. We have also made predictions for 
multistrange baryon production in heavy ion collisions at RHIC and found a similar enhancement. 
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I. INTRODUCTION 



One possible signal for the quark-gluon plasma (QGP) 
that is expected to be formed in relativistic heavy ion 
collisions is enhanced production of strange particles 
particularly those consisting of multistrange quarks such 
as cascade and omega as well as their antiparticles. The 
argument is that the rate for strange hadron production 
is small in hadronic matter due to the large threshold 
and small cross sections, which is contrary to the large 
production rate for strange quarks in the quark-gluon 
plasma For the lightest strange hadrons such as 

kaons and antikaons, experiments have shown enhanced 
production in heavy ion collisions at all energies avail- 
able from SIS (~ lA GeV) [|, AGS (~ lOA GeV) [|, 
and SPS 158A GeV) [|||. Studies based on transport 
models show that this enhancement can be explained 
by hadronic scattering alone However, these mod- 
els have failed in accounting for the enhanced produc- 
tion of multistrange baryons (S, fl) by about a factor of 
three, and the discrepancy is even larger for their antipar- 
ticles 0. To increase the production of these particles 
requires either to increase the string tension or to de- 
crease the constituent quark mass in the fragmentation 
of the initial strings in the dense matter |^] . It has also 
been suggested that the observed enhancement of mul- 
tistrange baryons and antibaryons may be due to topo- 
logical defects arising from the formation of disoriented 
chiral condensate in the initial high density stage of col- 
lisions Of course, the enhancement could simply be a 
result of the formation of the quark-gluon plasma during 
collisions. However, to establish strangeness enhance- 
ment as a signal for the quark-gluon plasma, we need 
to exclude any conventional mechanisms. As recently 
suggested by Vance strangeness-exchange reactions 
between antikaons and hyperons as well as between their 
antiparticles can contribute significantly to the produc- 
tion of multistrange baryons and antibaryons in relativis- 
tic heavy ion collisions. 

In this paper, we shall determine the contribution to 
multistrange baryon production from the hadronic mat- 
ter in heavy ion collisions at both SPS and RHIC en- 
ergies, using a recently developed multiphase transport 
(AMPT) model 0. In Section ||, we describe briefly the 



multiphase transport model. The strangeness-exchange 
reactions are then discussed in Section HI . Results from 



the AMPT model are given in Section IV B for heavy ion 
collisions at SPS and in Section IV C for heavy ion col- 



lisions at RHIC. Finally, a summary is given in Section 



II. MULTIPHASE TRANSPORT MODEL 



In the AMPT model, the initial conditions are obtained 
from the HIJING model ||ll[] by using a Woods-Saxon ra- 
dial shape for the colliding nuclei and including the nu- 
clear shadowing effect on minijet partons via the gluon 
recombination mechanism of Mueller- Qiu |l2|. After the 
colliding nuclei pass through each other, the Gyulassy- 
Wang model is then used to generate the initial 
space-time information of partons. Subsequent time evo- 
lution of the parton phase-space distribution is modeled 
by Zhang's Parton Cascade (ZPC) which at present 
includes only the gluon elastic scattering. After minijet 
partons stop interacting, they combine with their par- 
ent strings and are then converted to hadrons using the 
Lund string fragmentation model [|l^,^ after an aver- 
age proper formation time of 0.7 fm/c. Dynamics of the 
resulting hadronic matter is described by a relativistic 
transport model (ART) 0. 

The parameters in the AMPT model are fixed using 
the experimental data from central Pb-|-Pb collisions at 
center of mass energy of 17 A GeV Specifically, we 
have included in the Lund string fragmentation model the 
popcorn mechanism for baryon-antibaryon production in 
order to describe the measured net baryon rapidity distri- 
bution. Also, to account for the pion and enhanced kaon 
yields, we have modified the values of the parameters in 
the string fragmentation function. The same parame- 
ters are then used to study heavy ion collisions at RHIC 
energies. The model is able to describe the PHOBOS 
data ||l^ for the pseudorapidity distribution of charged 
particles Also, the predicted p/p and K~/Tr~ ra- 

tios are consistent with the data from the STAR [Ell and 



BRAHMS 1 22 1 Collaborations. Furthermore, it gives an 
elliptic flow |p3| which is comparable to that measured 
in the STAR experiment H]. 
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III. STRANGENESS-EXCHANGE REACTIONS 

To include multistrange baryon and antibaryon pro- 
duction, we consider the following reactions: 



KK .nTr, KYi <-> ii7r, and fiTr. 



(1) 
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FIG. 1. Experimental cross sections for K~p 
K~p E^tt" and K~n — > E°7r~ as a function of center 
of mass energy of the interacting antikaon and nucleon. 

Since there is no experimental information on these 
cross sections, we assume that they are the same as the 
cross section for KN Stt, which is the isospin averaged 
cross section for converting a nucleon to a sigma and can 
be related to the cross sections for K^p — > E^tt^, and 
K-n YPtt- by 



(2) 



The cross sections on the right hand side are known 
empirically and have been parameterized in Ref. [ p5| as 

follows. 



•0^0 = 0.6p-^-^ [mb]; 0.2 < p < 1.5 GeV/c, (3) 



1.2p-i-3 [mb], 

if 0.5 < p < 1 GeV/c; , , 

if 1 < p < 6 GeV/c. 



In the above, p is the K~ momentum in the laboratory 
frame. These cross sections are shown in Fig. |] as func- 
tions of center of mass energy of the interacting antikaon 
and nucleon. 

Empirically, the cross section for K^p Att" is also 
known and has been parameterized [M as 



50p2-67p-f 24 [mb], 

if 0.2 <p < 0.9 GeV/c; 

if 0.9 < j3 < 10 GeV/c, 

As shown in Fig. |l|, cri<--p^A^o has a similar magnitude 
as crx-p^T,°Ti° E^nd (JK-n^s^Tr- except that it has a res- 
onance contribution at ^/s ~ 1.8 GeV. It may be noted 
that the channel K^p Att^ has resonances while 
has no resonance and the branching ratios of the decay 
channels for the resonances of !B are not well determined. 
Therefore for a conservative estimate we assume that the 
reactions in Eq. (]l|) are similar to KN — > Stt, which does 
not show any resonance structure. 

We believe the assumption that all strangeness- 
exchange reaction cross sections have similar magnitudes 
is reasonable as these reactions are similar except for dif- 
ferences in their particle masses. Justification of this as- 
sumption is being studied in an effective hadronic model 
p6[ that is based on SU(3) flavor symmetry. Prelimi- 
nary results indeed show that the cross sections for these 
reactions have comparable values. 



IV. MULTISTRANGE BARYON PRODUCTION 
FROM HEAVY ION COLLISIONS 

A. Perturbative approach 

Since the abundance of multistrange baryons is small 
in heavy ion collisions, they can be treated perturbatively 
in the AMPT model as in kaon production from heavy 
ion collisions at low energies , where the collision dy- 
namics is assumed not to be affected by the production 
of these particles. 

First, we consider the process KA — > Stt, in which both 
K and A are treated explicitly, i.e., nonperturbatively, in 
the AMPT model. The KA collision is treated by taking 
their total scattering cross section to be U/fA = 20 mb. 
K and A then make a collision if their impact parameter 
is less than \fo^^jJ^- When this occurs, a S is produced 
with a probability given by Ph = R k-^b.-k I R Ki which is 
stored as the S formation probability. At the same time, 
both K and A have a probability of Ph being destroyed, 
and a pion is also produced with the same probability. 
Otherwise, the momenta of K and A are changed accord- 
ing to an isotropic distribution. 

To treat the inverse reaction of S absorption by a pion, 
i.e, Stt ^ KA, we use the cross section determined by de- 
tailed balance from Eq. (||). If a collision occurs, the S is 
destroyed. The pion has a destruction probability given 
by the S formation probability P-= and simultaneously 
a K and A are produced with an isotropic momentum 
distribution. 

The above perturbative approach can be similarly used 
for the reaction KH <-!■ Stt. However, it needs to be mod- 
ified for the reaction KB. ^ ilvr as both S and are 
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treated perturbatively. Again, we take the KE scatter- 
ing cross section to be cr^f^; = 20 mb. Whenever they 
colhde, the S have a destruction probabihty given by 
<7 Ks-*nTr / KS while an is produced with a probabil- 
ity given by the above probability multiplied by the S 
formation probability Ph, which is then stored as the f2 
formation probability Pn . The destruction of K and si- 
multaneous production of tt has a probability of Pq . The 
inverse reaction fin KE can be similarly treated as for 
En KA. 



proach, it was demonstrated |g8| that multimesonic re- 
actions Y + N WK -f nyK may lead to an enhanced 
antistrange hyperon Y production in a purely hadronic 
scenario. In the present transport model we have ne- 
glected these multiparticle reactions which in conjunction 
with the strangeness-exchange reactions could lead to a 
better reproduction of the experimental data for mul- 
tistrange baryons. Alternate to the multiparticle reac- 
tion, as for example, Stt -t- 2K ^ S + iV, the process 
k{k*) + K{k*) ^ E + N may be used instead in the 
transport model to generate S. 



B. Multistrange baryon production at SPS 
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FIG. 2. Time evolution of midrapidity hadrons for Pb+Pb 
collisions at SPS energy of ^/s = 17 A GeV at an impact pa- 
rameter of 6 < 3 fm in the AMPT model. 

In Fig. |[ we show the time evolution of the abun- 
dance of midrapidity kaons (including K*), antikaons (in- 
cluding K*), A, E, S^, and in Pb-|-Pb coUisions at 
^/s — 17 A GeV at an impact parameter 5 = — 3 fm. As 
is evident from the figure, most multistrange baryons are 
produced within 10 fm/c after the initial contact of the 
colliding nuclei when the energy density is high. Com- 
pared to the initial yield obtained from HIJING, the dy- 
namical evolution of the system leads to about 30% and 
70% increase in the A and production, while most of 
the Q^s are produced from hadronic rescattering. There- 
fore the observed enhancement of multistrange baryons 
can be largely accounted by the strangeness-exchange 
reactions in the hadronic matter. Furthermore, the fi- 
nal yield of strange baryons in the AMPT model are in 
reasonable agreement with that obtained by the WA97 
Collaboration The fl~ production in our model is 
however a factor of two smaller compared to the data. 
It may be mentioned that using a rate equation ap- 
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FIG. 3. Production and absorption rates of strange parti- 
cles in heavy ion collisions at SPS as functions of time. 

To see if strange hadrons reach chemical equilibrium 
in the collisions, we show in Fig. ^ the time evolu- 
tion of their production (solid lines) and absorption rates 
(dashed lines) . It is seen that most strange particles such 
as the hyperons and especially the kaons and antikaons 
approach equilibrium as their production and absorption 
rates become comparable before the system freezes out. 
So the abundance of these particles in heavy ion collisions 
at SPS are consistent with that given by the equilibrium 
models 12? 



In Fig. M we depict by open circles the transverse mass 
spectra of hadrons at midrapidity obtained in the AMPT 
model for the SPS energy. These results are obtained by 
using mean-field potentials for strange baryons that are 



3 



related to the nucleon potential according to their light 
quark contents. Specifically, relative to the nucleon po- 
tential, which is taken to be a stiff Skyrme potential with 
compressibility of 380 MeV, the lambda and sigma poten- 
tials are 2/3, the cascade potential is 1/3, and the omega 
potential is zero. We see that the theoretical transverse 
mass spectra of the hadrons agree well with the exper- 
imental data shown by solid circles |31 3^. The kaon 
transverse mass spectrum is however somewhat steeper 
than the experimental one. This may be due to the weak 
kaon mean-field potential at high density that is based on 
the impulse approximation Jl7t and is used in the present 
study. 
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FIG. 4. Transverse mass spectra for midrapidity protons, 
kaons, antikaons, pions (left panel), and lambda, cascade, 
omega (right panel) in 158^1 GeV Pb-|-Pb central collisions. 
The solid circles are the experimental data for ~ 5% most cen- 
tral collisions from the NA44 Collaboration [31] (left panel), 
and ~ 4% central events from the WA97 Collaboration [32] 
(right panel). The open circles are the AMPT model calcula- 
tions for impact parameter of 6 < 3 fm. 



The transverse mass spectra can be approximately fit- 
ted by an exponential function exp{—mT/T), where T is 
the inverse slope parameter. The slope parameter may 
provide useful information on the thermal motion and 
the strong collective transverse flow observed at SPS en- 
ergy. In Fig. 1^, we show the inverse slope parameter 
of these particles as calculated from the AMPT model 
(open circles) and compare them with those extracted 
from the experimental data (solid circles) The agree- 
ment between theory and experiment is reasonable. The 
linear increase of the inverse slope parameter from pion 
to kaon and to proton is indicative of the development of 
collective transverse flow in the system. For nonstrange 
hadrons, T may be parametrized by T = Tfo + m{/3t)'^ 
(dotted line), where the freeze-out temperature Tfo « 145 
MeV and average transverse flow {Pt) — 0.39c. As in 



the experimental data, the inverse slope parameters for 
strange baryons in the AMPT model are smaller than 
that for the proton. This is due to both the fact that 
strange hadrons are produced earlier in the collision when 
the energy density is high and have smaller scattering 
cross sections with pions and nucleons than nonstrange 
hadrons and thus freeze out quite early . Note that for 
proton with slope parameter T = 0.29 GeV, the contri- 
bution from HIJING is about 0.16 GeV, while the mean- 
field causes a further increase by « 0.03 GeV. 
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FIG. 5. Inverse slope parameters of various particles in 
heavy ion collisions at SPS. The open circles are the AMPT 
calculations and the solid circles are the experimental data 
[32] . The dotted line represents a linear parameterization (see 
text). 



C. Multistrange baryon production at RHIC 

The predictions of the AMPT model for multistrange 
baryon production in Au-t-Au collisions at RHIC energy 
of y/s — 130A GeV are given in Fig. ^ for their time 
evolution. Compared to that from heavy ion collisions at 
SPS, we see that the abundance of the multistrange par- 
ticles at this RHIC energy is increased by about a factor 
of two. We have also examined the time evolution of the 
production and absorption rates of strange particles. It is 
found that both kaons and antikaons approach chemical 
equilibrium. On the other hand, strange baryons such as 
S and 2 are somewhat out of equilibrium as their pro- 
duction rates are larger than their absorption rates at 
freeze out. Again, including strange baryon production 
from multiparticle reactions would lead to chemical equi- 
librium for these particles as well. Because of its larger 
absorption cross section by the pion compared to other 
strange baryons, the does reach equilibrium during col- 
lision. 
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We note that our predicted ratios of multistrange 
baryons to mesons, such as /K~ and r2~/7r~, are sim- 
ilar to that predicted by the thermal model |Q . In such 
statistical model analysis of the experimental data at the 
RHIC energy of -^/s — V6QA GeV, the hadron yields are in 
thermal and chemical equilibrium with a baryon chemical 
potential of /is ~ 51 MeV and freeze-out temperature of 
T ~ 175 MeV. The strangeness-exchange reactions intro- 
duced in the present study are thus another mechanism 
for multistrange baryons to achieve chemical equilibrium 
in relativistic heavy ion collisions. 
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FIG. 6. Time evolution of midrapidity hadrons for Au-|-Au 
collisions at RHIC energy of ^/s = 130A GeV at an impact 
parameter of 6 < 3 fm in the AMPT model. 
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FIG. 7. Same as Fig. ^ but for Au-|-Au collisions at 
^/s — 130A GeV for impact parameter of 6 < 3 fm in the 
AMPT model. 



In Fig. 1^ we show the transverse mass spectra of midra- 
pidity hadrons for central Au-|-Au collision at the RHIC 
energy y/s = 130A GeV. The spectra for protons, pi- 
ous and kaons are consistent with the preliminary data 
from the PHENIX Collaboration The agreement 

with this data is clearly evident in Fig. ^ where the 
slope parameters obtained by fitting the spectra in the 
same pt range with the exponential function is depicted. 
The mass dependence of T is qualitatively similar to that 
at SPS energy. At the RHIC energy the slope parame- 
ter from pion to proton also exhibits a strong depen- 
dence on the particle mass with the freeze-out tempera- 
ture Tfo ~ 155 MeV and the average collective velocity 
{j3t) ~ 0.43c. As at SPS the slope parameter for strange 
and multistrange baryons reveals a plateau since these 
particles, mostly generated by strangeness-exchange re- 
actions, are weakly interacting and decouples rather early 
from the system. 
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FIG. 8. Inverse slope parameters of various particles in 
heavy ion collisions at RHIC. The open circles are the AMPT 
calculations and the solid circles are the experimental data 
[35] . The dotted line represents a linear parameterization (see 
text). 



V. SUMMARY 



In this paper, we have used a multiphase transport 
model, that includes both an initial partonic matter and 
a final hadronic matter, to study the production of mul- 
tistrange baryons from the strangeness-exchange reac- 
tions in the hadronic matter. The cross sections for 
these reactions are assumed to be the same as the em- 
pirically known cross sections for antikaon-nucleon to 
hyperon-pion. For heavy ion collisions at SPS energies, 
we find that these reactions lead to an enhanced pro- 
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duction of multistrange particles, comparable to that ob- 
served by the WA97 collaboration. For heavy ion colli- 
sions at RHIC, a similar enhancement is obtained from 
om' model. We further find that the slope parameters for 
the (multi-)strange baryons exhibit a plateau as a func- 
tion of particle mass at both SPS and RHIC energies. 

We have not considered in the present study the 
hadronic contribution to multistrange antibaryons. For 
such a study, we need to extend the AMPT model to 
treat antibaryons similar to that for baryons, which is 
currently under way. 
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